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1. Introduction 

Strange quarks provide a special perspective on nucleon structure due to their absence from 
the valence content of the nucleon. As a result of this absence, effects of quark-antiquark pair 
creation and annihilation can be studied in isolation. Correspondingly, considerable efforts to probe 
strange quarks in the nucleon have been, and continue to be made on the experimental side. On the 
one hand, electron scattering experiments have provided data on the strange electric and magnetic 
form factors |JIJ], as well as given some indication of the strange axial structure; the latter, on the 
other hand, can also be studied in more detail in neutrino scattering experiments, with a combined 
analysis having been presented in [Q]. Considerably enhanced data are expected to emerge from 
the upcoming neutrino experiments MicroBooNE and MINER vA. 

On the side of lattice QCD, calculating strange contributions to nucleon structure counts 
among the relatively hard problems, since these contributions are determined purely by discon- 
nected diagrams, requiring propagator traces, cf. Fig. [l[ On the other hand, among this hard class 
of problems, strange matrix elements are still the easiest to access, since the quark being propa- 
gated in the disconnected loop is heavy, and the associated propagator calculation is therefore less 
expensive than for light quark loops. Thus, besides the physical relevance of such calculations, 
they can also serve as initial test cases for exploring techniques to evaluate disconnected contribu- 
tions to hadron physics more generally. A number of such investigations are being pursued by a 
variety of groups [j|-[TT]]. The effort presented here focuses on the two most fundamental contribu- 
tions of strange quarks to the properties of the nucleon, namely, their contributions to the nucleon 
mass and the nucleon spin. The calculational scheme employed is one which has been developed 



and optimized by the LHP Collaboration [|12|]. Domain- wall quark propagators are evaluated on 
HYP-smeared MILC asqtad dynamical quark ensembles. While the use of domain wall fermions 
implies considerable computational expense, it is expected to yield advantages in terms of mild 
renormalization and chiral behavior. 



2. Strange matrix elements 



The strange contributions to nucleon mass and spin can be characterized by the matrix ele- 
ments 

117 

f Ts = —(N\ss\N) and As = {N,i\s ri y 5 s\N,i) (2.1) 

m N 

respectively, where \N,i) denotes a nucleon state with spin polarized in the i direction. These 
matrix elements are obtained from corresponding lattice correlator ratios, 



R[V HK ,r obl ](z,T) 



I™ c Ejf Np (x, T)N a (0, 0)1 • s p (y, x)s y (y, z) 



^^^(1,1)^(0,0) 



(2.2) 



where N,N denote (smeared) nucleon sources and sinks, the sums over spatial position x project 
onto zero momentum nucleon states, the standard minus sign accompanying the quark loop has 
already become explicit through the reordering of the strange quark fields, and the T matrices im- 
plement nucleon polarization and operator insertion structure. Specifically, the unpolarized nucleon 
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r,t = r 
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N,t = J 



N, t = 



Figure 1: Disconnected contribution to nucleon matrix elements. The nucleon propagates between a source 
at t = and a sink at t = T; the insertion of F = T" bs occurs at an intermediate time t = T. 



two-point function in the denominator of (2.2) is achieved using Y unpo1 = (1 + Ya)/2; on the other 
hand, the three-point function in the numerator results from evaluating the correlation between a 
nucleon propagator and a strange quark loop, as also displayed diagrammatically in Fig. |j. To 
obtain fj , one calculates 



m N 



(R [T unpol A]iTj) _ [VEV] ) = R{fTs} _^ ^ 



(2.3) 



in the limit T T S> 0, where, as indicated, the vacuum expectation value of the quark loop, 
[V£V] = (— Zjt Sy{y,x)sy{y,x)) is subtracted; the matrix elements (2A) are meant to characterize 
the strange content of the nucleon relative to the vacuum, which itself contains a strange scalar 
condensate. The nucleon mass can be extracted from the nucleon two-point function as a by- 
product of the calculation. On the other hand, to obtain As, one calculates 



-i.2-R[(-iYiY 5 /2)r mpol ,YiY5K*,T) = R{^} 



As 



(2.4) 



in the limit T S> z » 0. In this case, the corresponding vacuum expectation value vanishes, but it 
can nevertheless be calculationally advantageous to subtract this numerical zero in order to reduce 



statistical fluctuations. Note that, in (2.4), an average has already been taken over the expectation 
values obtained using nucleons polarized in the positive and the negative /-directions, respectively; 
averaging the corresponding projectors (1 =F /^7 5 )/2 (with a relative minus sign) leads to the first 
argument of R in ( |2.4[ ). In the numerical calculation, also the polarization axis i will be averaged 
over the three spatial directions in order to further improve statistics, cf. the description further 



below. The prefactor 2 in < \2.4\ ) is a normalization factor compensating for the fact that the un- 



polarized nucleon two-point function is always used in the denominator of the ratio (|2.2|), even 
when polarized nucleon states are used in the numerator. Finally, the prefactor (— i) incorporates 
the Wick rotation back to Minkowski space; it compensates for the additional factor i arising in 



the 75 matrix contained in the second argument of R in ( |2.4D when casting the calculation on the 
Euclidean lattice. 



3. Lattice setup 

Numerical work was carried out on two 2+ 1-flavor asqtad dynamical quark ensembles pro- 
vided by the MILC collaboration, corresponding to the pion masses m % = 356MeV and m % = 
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Figure 2: Setup of the lattice calculation. The nucleon source is located at lattice time t = 0. The operator 
insertion time T is averaged over t = 3, . . . ,7; accordingly, stochastic sources are distributed over the bulk of 
the lattice in this entire time range. The temporal position T of the nucleon sink is variable. 



495 MeV. These ensembles contained 448 and 486 20 3 x 64 lattices, respectively, with lattice spac- 
ing a = 0.124fm. The configurations were HYP-smeared for the purpose of this calculation. The 
nucleon two-point functions and strange quark loops were evaluated using domain wall quarks. The 
lattice setup employed is depicted in Fig. ||. To enhance statistics, the operator insertion time X in 
( 2.2) was averaged over five time slices, t = 3, . . . ,7 (where t = corresponds to the nucleon source 
position). To implement this average, complex Z(2) stochastic sources, introduced to the evaluate 
the strange quark loop propagator trace, were distributed over the bulk of the lattice within this en- 
tire temporal range. For m n = 356MeV, 1200 stochastic sources per configuration were used, for 
m % = 495 MeV, 600 stochastic sources 1 . This rigid temporal setup, chosen to allow for maximal 
statistics, does not permit a variation of the operator insertion time in order to test for a plateau; it 
is motivated by previous extensive experience with connected diagrams in the same scheme Jl2|], 
which suggests that the nucleon ground state has been filtered out at t = 3 to a sufficient degree as 
to render the associated systematic uncertainty small compared to the statistical uncertainty of the 
present calculation. On the other hand, a residual opportunity to test the dependence of the results 
specifically on the operator-sink separation is given; below, results for the relevant correlator ratios 
will be shown as a function of variable sink position T, with the expectation that the asymptotic 
behavior will be reached for sink positions T > 10. 

Besides statistical uncertainties stemming from the stochastic evaluation of the quark loop, the 
observables studied here exhibit substantial gauge fluctuations. To accumulate statistics sufficient 
to overcome these fluctuations, it is necessary to evaluate multiple samples of the correlator ratios 
of interest per given gauge configuration. Thus, for given source time slice, not one, but several 
(eight in the case of m K = 356 MeV, four in the case of m K = 495 MeV) different samples were 
obtained by varying the spatial source position. Furthermore, since the scheme depicted in Fig. || 
requires a much smaller temporal extent than available on the lattice employed, the entire scheme 
was replicated three times on separate temporal regions of the lattice. Altogether, therefore, 24 
correlator ratio samples were obtained per lattice gauge configuration in the case of m n = 356 MeV, 



'These high statistics in the stochastic sources mainly serve to improve the signal for As; the stochastic estimator 
for the scalar matrix element, by contrast, converges rapidly and the statistical uncertainty in fx s is dominated by gauge 
fluctuations. 
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Figure 3: Correlator ratio R{fr s } as a function of sink time T, for the two pion masses considered. 
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Figure 4: Correlator ratio /JjAs} as a function of sink time T, for the two pion masses considered. 

and 12 in the case of niji — 495 MeV. Finally, as already indicated in the previous section, As was 
averaged over three separate polarization axes i, corresponding to the coordinate axes. 



4. Numerical results and conclusions 

Figs. H and |4] display, as a function of sink time T, the correlator ratios R{fT s } and /?{As}, 
cf. (2.3) and ( |2.4| ), averaged in the fashion described in the preceding section; for large T, these 
quantities yield the strange contributions to the nucleon mass and spin, fj s and As. Before reaching 
the physically most relevant region T > T (where z denotes the temporal location(s) of the operator 
insertion), the ratios start out at vanishing values near T = 0, then gather up magnitude as the 
source-sink time interval includes an increasing portion of the stochastic source region contributing 
to the quark loop. In the region T > X, the correlator ratios are expected to level off to approach their 
asymptotic value for large sink times. On the the other hand, in this region, statistical fluctuations 
become appreciable, and there is therefore only a short time window in which one can hope to 
observe this behavior. 

In the case of As, cf. Fig. Q this expectation is confirmed rather well. Quantitatively, the 
estimates of As obtained either from the correlator ratio at sink time T = 10 or from the average 
over the correlator ratios at sink times T = 10, . . . , 14 yield almost identical results, 



At w% = 356MeV : 
At m n =495 MeV : 



As| r= io = -0.030(9) 
As| r= i = -0.030(7) 



As|r=io...i4 = -0.030(19) (4.1) 
Aj|t=io...14 = -0.031(11) (4.2) 
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Figure 5: Pion mass dependence of the results for fj s and As. Filled circles represent measured data; open 
circle in left-hand panel shows chiral extrapolation of the data to the physical pion mass, cf. main text. 



where the error estimate in the sink-time averaged case is the jackknife error extracted by per- 
forming the sink time average configuration by configuration, i.e., the fact that correlator ratios at 
different sink times are not independent is taken into account. 

On the other hand, in the case of fj s , the behavior of the correlator ratios as a function of sink 
time is not as clear-cut. At m K = 356 MeV, the correlator ratio for T > 10 considerably overshoots 
the value at T = 10, whereas at m n = 495 MeV, the correlator ratio for T > 10 decreases again 
compared to its T = 10 value. Quantitatively, the comparison analogous to the one performed for 
As above yields 

Atm^ = 356MeV: f Ts \ T=W = 0.043(6) /rJr=io...i4 = 0.057(11) (4.3) 

Atm^=495MeV: /r,|r=io = 0.046(4) /rJr=io...i4 = 0.037(5) (4.4) 

The fact that the deviations from the expected plateau behavior occur in opposite directions for the 
two pion masses may be an indication that they are caused by statistical fluctuations; also the error 
estimates in (|4.3[),(fk4]) are still compatible with this possibility. In view of the fact that correlator 
ratios at different sink times are not independent of each other, the rather smooth behavior of the 
correlator ratio in the case of m n = 495 MeV does not necessarily contradict an explanation in terms 
of statistical fluctuations. Currently, a doubling of the statistics for the m % = 495 MeV case is being 
pursued to further explore this issue. 

Fig. H summarizes the results obtained as a function of pion mass, using the correlator ratios 
at sink time T = 10 to estimate the observables; for the case of fj s , the above discussion of the un- 
certainties involved in identifying an asymptotic value of the corresponding correlator ratio should 
be kept in mind. For fj s , also a tentative extrapolation to the physical pion mass is displayed; 
since the strange scalar matrix element in the nucleon is related via the Feynman-Hellman theorem 
to the derivative of the nucleon mass with respect to the strange quark mass, the chiral behavior 
follows from differentiation of nucleon mass formulae obtained in chiral effective theory [|7|]. The 
leading dependence is linear in the light quark mass, i.e., in m 2 n ; this was used in the extrapolation 
shown in the left-hand panel of Fig. ||. The extrapolated value is fj s = 0.041(12), corresponding to 
m s (N\ss\N) = 39(12) MeV. 

On the other hand, while no attempt at extrapolating As in the right-hand panel of Fig. || has 
been made, no significant variation with pion mass is seen. It should be noted that the strange axial 
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current entering the calculation of As requires renormalization, which has not been included in the 
results presented here; however, for the lattice scheme used in this work, axial current renormal- 
ization constants translating to the MS scheme at a scale of 2 GeV are consistently very close to 



1.1, over a wide range of quark masses [J12p. It is therefore expected that the results for As only 
acquire a mild 10% enhancement when translated to the standard MS scheme. Thus, no evidence 
for unnaturally large strange quark contributions to nucleon spin is seen in the present calculation. 

Of course, it should be noted that no attempt has been made at this point to quantify several 
other sources of systematic uncertainty, such as the dependence on the lattice spacing, lattice size, 
and, in particular, the fact that the strange quark mass in the gauge ensembles used lies appreciably 
above the physical strange quark mass. Data on the strange quark mass dependence reported in 
[pH, translated to the present scheme, indicate that matrix elements such as (N\ss\N) acquire a 
correction amounting to about 15% (which in fj s is (over)compensated by the m s prefactor); it thus 
appears reasonable to conjecture a 15% systematic error also for As from this source, about half 
the magnitude of the statistical uncertainty. 
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